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Abstract. It is shown that every scalar linear quadrilateral lattice equation lies 
within a family of similar equations, members of which are compatible between 
one another on a higher dimensional lattice. There turn out to be two such 
families, a natural parametrisation is given for each. 



PACS numbers: 02.30.Ik 

1. Introduction 

By integrability we mean structure in non-linear systems, however it is quite natural 
for a specific criteria for integrability to also admit linear examples. From the point 
of view of classification it is desirable to exhibit these systems. There is also the 
possibility that the linear examples provide a useful context for methods of the non- 
linear theory. Or that the non-linear theory can shed light on linear problems, a point 
of view which has been successfully exploited for linear PDEs by Fokas and co-authors, 
for example see [1]. 

Loosely speaking, the multidimensional consistency of an equation means 
identifying the equation as lying within a parametrised family of similar equations, 
members of which are compatible between one another on a higher dimensional lattice. 
This property, which was first made explicit in [2] and later also in [3], is sufficient 
for the integrability of scalar quadrilateral lattice equations and has led, through 
pioneering work of Adler Bobenko and Suris, to classification results in this area 
[4, 5, 6]. In these works the classification problem is tackled through the associated 
classification of polynomials of degree one in four variables which are compatible 
when associated to the faces of a cube. This approach to classification also picks 
out compatible systems which are not symmetric between the various faces of the 
cube [6] (in [4, 5] additional symmetry was imposed and the non-symmetric cases 
were therefore excluded, non-symmetric systems on the cube were further discussed 
in [7, 8]). 

However, a remarkably strong result was formulated and proved in [6] by 
introducing a notion of non-degeneracy for the defining polynomial, such polynomials 
were termed 'type-Q'. It was shown that every such polynomial lies in a unique 
(parametrised) compatible family of similar type-Q polynomials. The classification 
question for scalar quadrilateral lattice equations based on the multidimensional 
consistency is not fully answered for equations defined by polynomials which are not of 
type-Q. Certainly there are non-trivial examples, the polynomial defining the lattice 
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potential KdV equation [9, 10] is itself not of type-Q. All linear equations also lie 
outside this class and it is these equations that are studied in the present paper. 

We proceed by asking for compatibility between generic autonomous scalar linear 
quadrilateral lattice equations on a three dimensional lattice, this yields an algebraic 
system which constrains the coefficients of the equations. It turns out that this system 
is tractable by methods from linear algebra. Using the appropriate transformation 
group and a re-parametrisation we then suggest a canonical form for the compatible 
systems which emerge. 

2. Imposing consistency on the cube 

We will ask for compatibihty between the following system of polynomial equations, 

aiu + b\u + ciu + u = di, aiu + biu + ciu + u = di, 

a2U + b2U + C2U + U = d2, a2U + b2U + C2U + U = d2, (1) 

asu + bsu + C3U + u = c?3, aau + b^u + c^u + u = ^3. 

Here ai,bi,Ci 6 C \ {0} and G C are the coefficients of the polynomials. By 
compatibility we mean that given initial data u, u, u and u taken from C, and after 
evaluation of the intermediate variables u, u and u using the equations on the left, 

the remaining equations on the right determine the same value for u. By assigning 
the variables involved to the vertices of a cube as in Figure 1 each equation in (1) 
may be associated to a face of the same cube. This provides a convenient geometrical 
configuration to visualise this notion of consistency as well as giving the property its 
name [2, 3]. The system of equations (1) is not the most general system of linear 




u 



Figure 1. Variables associated to the vertices of a cube 

polynomial equations on the cube because equations associated to opposite faces 
coincide. However, this restriction is natural because polynomial equations which are 
compatible in the sense described may be used to define an autonomous compatible 
three-dimensional lattice system (as in the next section). 

It is a matter of straightforward calculation to verify that the system (1) is 
compatible as described if and only if the coefficients satisfy the following system 
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of equations: 

bibs - 63C2 + C2C1 = ai, 

62&I - ^lC3 + C3C2 = 02, (2) 
&3&2 — + C1C3 = 0,3, 

(62-C3)di - (1 + 63)^2 + (1 + C2)d3 ^ 0, 

(1 + C3)di + (fe3-ci)d2 - (1 + 6i)d3 = 0, (3) 

-(1 + 62)^1 + (1 + Cl)d2 + (fol-C2)d3 0. 

In tlie remainder of this section we give a prehminary analysis of this system. In 
particular we will separate out three cases based on the rank of the matrix emerging 
in (3) which we denote by M, 

62 - C3 -1 - &3 1 + C2 

M := I + C3 63 -ci -I-&1 . (4) 

-1 - 62 1 + Ci 61 - C2 

There are only three cases because rank(M) < 2, so for example the last equation 
in (3) is redundant because it is a consequence of the first two. Suppose first that 
the rank of M is zero. This happens only if all entries of M vanish, from which the 
following is immediate, 

rank(M) = <^ 6^ = q = -1, iG {1, 2, 3}. (5) 

Consider now the case that rank(M) = 1. If ai, 02 and 03 are fixed in terms of bi, Ci, 

i e {1,2,3} by (2), then 

rank(M) < 1 <^ a, + h + c, + 1^0, {1,2,3}. (6) 

This is true because the condition + + Cj + 1 = 0, i G {1, 2, 3} is equivalent to 
the vanishing of all cofactors of M which is necessary and sufficient for rank(M) < 1. 
The following is also true, 

rank(M) = 1 kcrncl(M) = span(77i, 112, U3), (7) 

where the vectors vi, V2 and U3 are defined as 



(8) 



Assuming the condition on the right of (6) holds it can be verified by calculation that 
the vectors V2 and W3 defined in (8) lie in the kernel of AI . and by further calculation 
that span((,'i. '(;2! V3) has dimension less than 2 only if hi = Ci — —1, i G {1,2,3}. 
But by (5) this latter condition violates rank(M) = 1, hence if rank(M) ~ 1 then 
span(t)i, f2, U3) must have dimension 2, which completes the verification of (7). And 
finally, 

rank(M) = 2 kernel(M) = span(«;), (9) 

where the vector w is defined as 

ai +bi + ci + I 

w:= 02 + 62 + C2 + 1 , (10) 
_ as + ^3 + C3 + 1 _ 

and again ai, 02 and 03 are given in terms of bi, Ci, i G {1, 2, 3} by (2). A calculation 
shows that w lies in the kernel of M, whilst by (6) the vanishing of w contradicts 
rank(M) = 2 so span(u;) has dimension 1 and (9) is verified. 
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3. Suggested canonical forms 

We consider here the three dimensional lattice system 
aiu + biu + ciu + u = di, 

a2U + b2U + C2U + U = d2, (11) 
asu + 63U + C3U + u = ds, 

which is similar to the polynomial system (1), but where now u ~ u{n,m,l), 
u = u{n + 1, m, I), u = u{n,m + 1, 1), u = u{n, m, i + 1) and u = u{n + 1, m + 1, 
etc. are values of a scalar dependent variable w as a function of three independent 
variables n,m,l £ Z. The coefficients a^, bi, Ci ^ C \ {0} and di E C, i E {1, 2, 3} are 
assumed to be autonomous. By the compatibility of (11) we will mean compatibility 
of the defining polynomials as described in the previous section, specifically that the 
coefficients satisfy the equations (2) and (3). 

It is clear that the form (and compatibility) of (11) is preserved by the point 
transformation group T 

T = { u{n, m, I) au{n, m, Z) + /3 | a e C \ {0}, /3 e C } . (12) 

And if the coefficients are such that (11) is invariant under the action of an abelian 
subgroup S<T, then the form of (11) is also preserved by the transformation group S, 

S = { u{n, m, I) [x" • • z'-]{u{n, m, 0) I x, y, 2 e S } , (13) 

elements of which are often referred to as gauge transformations. (If (11) were not 
invariant under S, then transformations from S would render it non- autonomous.) In 
the following theorem we use these transformations to suggest a canonical form for 

compatible systems (11). 

Theorem. By the action of point and gauge transformations, compatible systems (11) 
may be brought either to the form 

u — u — u + u = di, 

u — u — u + u = d2, (14) 
u — u — u + u = da, 

for some choice of the parameters di, ^2, c?3 € C, or to the form 

{pi - qi)u + (p2 - q2)u = {p2 - qi)u + (pi - q2)u, 

{qi -ri)u+ {q2 - r2)u = (52 - ^-i)" + (qi - r2)u, (15) 

(n -pi)u+ (r2 - P2)u = (r2 -pi)u+ (ri - p2)u, 

for some choice of the parameters Pi,qi,ri G C, i € {1,2}. Conversely, either of 
the systems (14) or (15) are compatible for any choice of their respective parameters 
di,d2,d3 e C and Pi, qi, e C, i G {1, 2}. 

Proof. To the system (11) we associate a matrix M defined by (4). Compatible 

systems (11) for which M is rank zero are necessarily of the form (14), a fact which is 
evident by inspection of (5), (2) and (3). We will now demonstrate that by point and 
gauge transformations, any compatible system (11) for which the associated matrix M 
is of non-zero rank can be brought to a form where rank(M) < 1 and di = ^2 = ^3 = 
(we refer to systems satisfying this latter condition as homogeneous). 
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Consider first that rank(M) = 1. Using (7) we know that di, d2 and ds are 
necessarily of the form 

di =/3i(l + 6i)+/32(l + ci), 

d2=/32(l + 62)+/33(l + C2), (16) 
d3=/33(l + 63)+/3l(l + C3), 

for some Pi, 02,^3 € C. But the condition on the coefficients appearing in (6) is 
equivalent to the invariancc of (11) imdcr transformations of the form u{n,m,l) — > 
u{n, m,l) + (3 for any /3 e C. It follows that we may apply the gauge transformation 

u{n, m, I) u{n, m, I) + + mf32 + Wa, (17) 

which brings (11) to the homogeneous form. Note that M does not change under this 
transformation so that in particular we still have that rank(M) = 1. 

Consider now the only other case, that rank(M) = 2. It is immediate from (9) 
that di , d2 and d^ are necessarily of the form 

di = S{ai + 6i + ci + 1), 

d2=5{a2+b2 + C2 + l), (18) 
^3 = 6{as + 63 + ^3 + 1), 

for some (5 S C. Applying the point transformation u{n, m,l) ^ u{n, m, I) +6 therefore 
brings (11) to the homogeneous form. This transformation does not change M, 
however it docs leave (11) with an invariance to point transformations of the form 
u{n, m, I) au{n, m, ^) for a e C \ {0}. We now claim that application of the gauge 
transformation 

u{n, m, I) aia^a^gu{n, m, I), (19) 

where the parameters ai, a2 and as are defined in terms of the coefficients and a new 
parameter fi;eC\{— l,l}by the equations 

-2 1-K, 
ai = — — ci + 63, 

1 + K 1 + K 

1 + K 2 

a2 = C2 - 61, (20) 

1 — K 1 — K 

K — 1 1 + K, 

"3 = C3 ^02, 

brings the homogeneous case of (11) to a form in which the associated matrix M is at 
most rank 1. The claim can be verified by calculation, the calculation is equivalent to 
verifying that homogeneous systems (11) admit the solution u{n,m,l) = a"a2*a3'Uo 
where uq G C is an arbitrary constant. This follows because homogeneous systems (11) 
admit non-zero constant solutions only if the condition on the right of (6) holds. Note 
that to ensure invertibility of the transformation (19) we should choose k G C\{ — 1,1} 
so that none of ai, 0:2 or 0:3 vanish, inspecting (20) we see this can always be done, 
pathological cases are avoided because bi, Ci, i G {1,2,3} are assumed non-zero. 

So up to point and gauge transformations, we may assume that compatible 
systems (11) are of the form (14) or else are homogeneous and satisfy the condition 
on the right of (6). It remains to show that systems of the latter form admit the 
parametrisation in (15). 
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This final step in the proof is again through calculation, choose distinct parameters 
P2,<12,1~2 G C and define pi, qi and ri by the equations 

P2 - 92 r2-p2 , 

Pi = T2C1 H ^203, 

92 -r2 92 - r2 

92 - r2 , P2 - 92 , 
9i = P2C2 H r26i, (21) 



P2 


- 92 


92 


- r2 


92 


- r2 


r2 


-P2 


r2 


-P2 


P2 


- 92 



r2 


-P2 


92 


- r2 


P2 


- 92 


r2 


-P2 


92 


- r-2 



ri = g2C3 H P2O2. 

P2 - 92 

Substitution of (21) into (15) followed by use of the condition on the right of (6) and the 
relations (2) brings (15) back to parametrisation in terms of the original coefficients. 
The second statement of the theorem can be verified by calculation. 

□ 

The compatible system (15) has been re-written in terms of lattice parameters, 
these are the two-component vectors (pi,P2), (91192) and (ri,r2). By fixing an 
association between these three parameters the three lattice directions n, m and I 
respectively the system has an important property. Namely covariance, it is invariant 
under permutations of the lattice directions. The natural extension of the compatible 
system to any higher number of dimensions becomes apparent from this property. 

We remark that the extension of (14) to higher dimensions is also clear, but in 
this case it is natural do associate the constants di each to a pair of lattice directions. 
This has the consequence that the number of constants present in the system grows 
quadratically with the dimension of the system as opposed to linearly in the case of 
(15). 

We make a historical remark regarding the parametrisation of the compatible 
linear system (15). The parametrisation (15) appeared in [7] in connection with a 
multidimensionally consistent equation discovca'od by Hictarinta [11]. The Hietarinta 
equation was shown to be lincarisable originally in [12] with a different parametrisation 
given for the associated linear equation. In the special case p2 = —pi, 92 — —Qi and 
r2 = — ri the system (15) reduces to the equation found by linear approximation of the 
NQC equation [13], which is the equation satisfied by the discrete soliton plane- wave 
factors denoted pk in [13]. Note that the corresponding solution of (15) is 

Finally we note that the two dimensional equation appearing in the homogeneous 
case of system (14) was studied as an example in [14]. 

4. Discussion 

It is easily verified that by the action of point and gauge transformations, quadrilateral 
lattice equations of the form 

au + bu + cu + u = d, (23) 

where a,b,c £ C \ {0}, d & C and u = u{n, m), u = u{n + 1, m) etc., may be brought 
to one of the two forms 

u-u-u + u = d, deC, (24) 

u-u = ^{u-u), 7 e C \ {0}. (25) 
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Comparing this observation with the systems (14) and (15) occuring in the main result 
reveals that all scalar linear quadrilateral lattice equations are multidimensionally 
consistent. But note that (24) has more potent solvability property beyond its linearity 
and multidimensional consistency, specifically for solutions of (24) it is true that 

u{0, 0) — u{n, 0) — u(0, m) + u{n, m) = nmd, (26) 

i.e., the equation can be directly integrated. The principal example of the linear 

multidimensionally consistent equation is perhaps then the equation (25), and on 
comparing this with equations in the system (15) there appears to be some redundancy 
in the three dimensional system. However, choosing for example 

Pi = 1-7, 91 = -1-7, P2 = -l + 7, 92 = 1+7, (27) 

in (15) removes the redundancy, although it breaks the covariance. Actually the 
covariance is still present, but hidden: Treating the system as an equation and its 
Backlund transformation (with Backlund parameter (ri,r2)) the covariant system 
emerges again as the superposition principle. 
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